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SUMMARY 


Weight  efficiency  when  using  thin  plates  in  aircraft  structures 
requires  increasing  load  carrying  capacity  beyond  buckling  load.  In 
this  respect  flat  plates  supported  along  their  edges  differ  favourably 
from  columns  and  curved  plates. 

The  aim  is  to  find  the  relation  between  the  edge  displacements  of  a 
panel  and  its  load,  including  the  relation  between  the  increments  of 
both  buckling  and  stiffness. 

The  physical  characteristic  feature  is  that  the  deflect ions  are 
finite,  from  which  the  equations  are  non-linear.  Even  with  simple 
boundary  conditions  exact  solutions  are  missing.  Approximate  solutions 
based  upon  energy  theorems  are  available.  Solutions  established  for 
flat  plates  at  constant  temperature  and  for  thermal  loading  will  be  dis¬ 
cussed  and  compared  with  available  experimental  evidence. 

Prom  a  recent  investigation  it  appears  that  the  post  buckling  behaviour 
of  narrow  cylindrical  panels  can  present  the  unstable  explosive  c-harac 
ter,  known  with  full  cylindrical  shells,  at  least  at  deformations  clos#.* 
to  those  at  which  buckling  starts. 
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SOMMAIRE 


Lj*  emploi  de  plaquettes  minces  demande  des  structures  capables  de  resist 
aux  charges  superieures  a  la  change  de  flambage  A  cet  egard  les  plaquettes 
minces  supportees  le  long  de  leurs  bords  presentent  des  caract^rtstiques 
favorables  vis-a-vis  celles  des  colonnes  et  des  plaquettes  courb6es. 

On  cherche  a  etablir  la  relation  entre  les  deplacementsdo  bord  d’ un 
panneau  de  revetement  et  la  charge  suhie,  y  enmpris  le  rapport  entre  line 
augmentation  de  flambage  et  de  rigidite. 

La  particularite'  pnysique  est  que  les  deflections  sont  finies,  donnant 
ainsi  des  equations  non  lineaires.  Meme  dans  des  conditions  de  limite 
simples,  des  solutions  precises  ne  sont  pas  possibles.  On  dispose  de 
solutions  fondees  sur  des  thcoremes  energe'tiques.  Des  solutions 
etablies  pour  des  plaquettes  planes  a  temperature  constanto  soumises  a 
des  sollicitations  thermiques  sont  etudi^es  et  comparees  avec  les  resul- 
fats  experimentaux  disponibles.  D'apres  une  etude  recente  il  parait 
quo  le  comportement  aprks  flambage  de  panneaux  cylindriques  etroits 
peuvent  presenter  le  caractere  instable  explosif  bien  connu  dans  le  ens 
de  coquilles  cylindriques,  tout  au  moins  lors  des  deformations  voisines 
k  celles  provoquant  le  flamhage. 


539.384.4 


2b2c3 


CONTENTS 


Page 


SUMMARY  i  i 

LIST  OF  FIGURES  v 

NOTATION  vi 

1.  INTRODUCTION  1 

2.  VARIETIES  OF  POST- BUCKLING  BEHAVIOUR  2 

3.  THE  PHYSICAL  CONCEPT  OF  POST  BUCKLING  BEHAVIOUR  4 

4.  THE  EQUATIONS  GOVERNING  POST-BUCKLING  BEHAVIOUR  OF 

FLAT  PLATES  7 

5.  METHODS  FOR  SOLUTION  9 

6.  AVAILABLE  SOLUTIONS  FOR  BEHAVIOUR  AFTER  BUCKLING  11 

6.1  Longitudinal  Compression  (Effective  Width)  11 

6.1.1  Flat  plates  11 

6.1.2  Narrow  cylindrical  panels  13 

6.2  Flnt  Plates  Under  Shear  Load  14 

6.3  Flnt  Plates  Under  Thermal  Stresses  17 

REFERENCES  20 

FIGURES  22 

DISTRIBUTION 


lY 


LI ST  OF  FIGURES 


Fig.l  Schematic  representation  of  equilibrium  load  versus  buckling 
deflection 

Fig. 2  Schematic  representation  of  equilibrium  load  versus  edge 
displacement 

Fig. 3  Schematic  representation  of  the  effect  of  imperfections 

Fig  4  Fcl>  ratic  representation  of  the  effect  of  initial  imperfections 

Fig. 5  Relative  decrease  of  membrane  stresses  with  increasing  edge 

displacement  beyond  buckling.  Longitudinal  compressive  strain  e 
of  a  long  rectangular  plate  of  width  b 

Fig. 6  Shortening  of  a  plate  due  to  buckling 

Fig. 7  Shear  angle  of  a  rectangular  plate  following  from  oblique  waves 

Fig. 8  Wave  pattern  for  large  strain  ratio  7/7^ 

Fig. 9  Stiffness  of  buckled  panels 

Fig. 10  Explosive  change  of  buckle  pattern 

Fig. 11  Effective  width 

Fig. 12  Effective  width  in  the  plastic  range  with  24S-T  (+)  find  75S-T  (0) 
cladded  and  uncladded  material 


Fig. 13  Effective  width  in  the  plastic  range  with  2S-&H  cladded  material 

Fig. 14  Tangent  to  post  buckling  curve  at  buckling  load  for  cylindrical 
panels 


Fig. 15  Effect  of  imperfections  on  maximum  load  in  the  initial  post 
buckling  range 


Fig. 16  Schematic  representation  of  post-buckling  behaviour  in  advanced 
stage  for  narrow  cylindrical  panels 


Fig. 17 
Fig. 18 


Direct  stresses  of  a  skin  plate  by  thermal  and  mechanical  load 


Temperature  distribution  in  the  plate  y  = 


Fig. 19  Load  strain  curve  with  thermal  conditions.  P,  e  applies  to  non- 
thermal  conditions.  P,  e  applies  to  thermal  conditions 


Page 


22 


22 


23 

24 


25 

26 
26 
27 
27 


28 


28 

29 

30 


31 


3) 

32 

33 
33 

33 


v 


NOTATION 


A 

i* 

R 

T 

b 

\ 

b 


11 

a 

y 

e 

0 


a 


y/7„ 


(aT)av 
P,  e 
q 


.“train  energy 
load 

radius  of  cylindrical  panel 
temperature 
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ns  a  suffix  denotes  the  quantity  concerned  at  buckling  load 

edge  displacement 

coefficient  of  thermal  expansion 

overall  shear  of  a  panel,  following  from  edge  displacement 

overall  longitudinal  compressive  strain,  edge  strain 

curvature  parameter  defined  by  Equation  (9) 

compressive  stress 
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P,  e  in  isothermal  condition 
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TOST  RUCKLING  HEHAVI«VJR  OF  STRUCTURES 


A.  van  der  Neut* 


1.  INTRODUCTION 

Engineers  have  for  a  long  time  been  Interested  In  the  buckling  of  structures,  but 
they  did  not  pay  attention  to  the  behaviour  of  structures  after  the  onset  of  buckling. 
The  accepted  Idea  was  chat  the  buckling  load  was  the  ultimate  performance  of  the 
structure,  that  for  reasons  of  safety  the  actual  load  had  to  be  kept  far  below  what 
was  called  the  critical  load  and  that  research  on  post-buckling  behaviour  had  no 
practical  importance. 

The  experience  that  some  structures,  such  as  axially  compressed  cylindrical  shells, 
failed  at  loads  far  below  the  theoretical  buckling  load  was  attributed  to  inevitable 
Imperfections  in  the  shape.  This  seemed  to  support  the  appreciation  of  critical 
loads  and  was  answered  by  the  Introduction  of  an  ample  factor  of  safety.  On  the 
other  hand  the  experience  that  other  structures,  for  instance,  plates  supported  at 
their  edges,  were  able  to  sustain  loads  far  in  excess  of  their  buckling  load,  was 
considered  with  suspicion  and  did  not  upset  the  belief  that  the  loads  should  be  kept 
below  the  critical  load. 

Aeronaut ic.J  engineers,  however,  could  not  close  their  eyes,  either  to  the 
treacherous  characteristics  of  shells  or  to  premising  phenomena  of  plates,  since  they 
were  living  under  the  ever  increasing  demands  of  weight  economy  in  structures  which 
were  highly  liable  to  buckling. 

When  applying  plates  as  structural  elements  of  aircraft,  these  plates  tend  to  he 
thin  and  consequently  buckling  stresses  are  much  smaller  than  the  allowable  strei  ses 
for  more  solid  sections  made  from  the  sane  material.  The  most  economical  solution 
is  to  rearrange  the  distribution  of  material  so  that  by  having  numerous  supporting 
elements  the  buckling  stress  of  the  reinforced  plate  is  well  up  to  the  material 
strength.  With  heavily  loaded  aircraft  structures  this  can  be  achieved,  as  in  the 
cases  of  the  main  structure  of  large  aspect  ratio  wings  and  of  thin  wings  for  large 
wing  loading  and  load  factor.  However,  with  plates,  where  the  load  per  unit  of  length 
is  relatively  small,  a  design  in  which  buckling  would  be  admitted  only  at  the  ultimate 
load  would  be  uneconomical. 

Such  conditions  are  present  with  lightly  loaded  wing  structures  and  the  majority 
of  fuselage  structures.  Here  plate  buckling  has  been  accepted,  in  many  cases  even 
in  horizontal  flight,  since  this  type  of  buckling  is  not  catastrophic  and,  moreover, 
the  plate  takes  its  share  in  supporting  load  Increments  due  to  manoeuvres  and  gusts. 
Such  a  plate  can  be  regarded  as  remaining  a  structural  element  even  after  buckling. 

To  exploit  this  favourable  characteristic,  aeronautical  engineers  hnd  to  investi¬ 
gate  this  ‘post-buckling  behaviour'  of  plates. 

No  knowledge  acquired  by  other  branches  of  engineering  was  available.  Although 
civil  engineering  was  confronted  with  the  same  type  of  prob1  ,  for  instance  deep 
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l-hcams  In  bridges,  the  need  for  weight  economy  was  not  so  pressing  there  that  post- 
hucklirtj;  use  was  acknowledged  to  p~csent  a  sound  engineering  philosophy. 

It  cannot  be  denied  that,  some  dawn  of  consciousness  of  post-buckling  theory 
existed  in  nuval  architecture.  Then  the  steel  skin  of  ships  was  reinforced  by 
longitudinals,  a  snail  width  of  plate  was  considered  to  he  effective  as  an  addition 
to  the  rection  of  the  longitudinal.  However,  due  to  the  Classification  Rules,  strength 
analysts  could  not  lead  the  way  to  weight  saving  structures  and  the  investigation  of 
post-buckling  behaviour  was  not  undertn'-en.  So,  as  in  other  cases,  to  aeronautics 
fell  the  honour  of  penetrating  the  brambl e-bashes  and  awnkening  Sleeping  Beauty  with 
its  scientific  kiss.  If  this  metaphorical  language  calls  for  a  persenif icatlon  there 
is  no  doubt  that  this  romantic  part  may  be  attributed  to  H.  Wagner,  who  published  in 
1929  his  wellknown  work  on  the  tension  field  of  spar  webs’.  This  paper  indeed  is  the 
charter  legitimating  the  use  of  buckled  plates  as  structural  elements  for  carrying 
shear  loads  in  spar  webs  and  in  the  skins  of  wings  and  fuselages. 

In  the  early  thirties  the  American  aircraft  industry  introduced  a  new  and  highly 
important  application  of  post-buckling  with  so-called  stressed-skin  structures. 
Supporting  the  skin  by  numerous  longitudinal  stringers  the  skin,  though  buckled,  could 
be  made  effective  in  carrying  longitudinal  compression.  Knowledge  of  the  behaviour 
of  this  type  of  structure  started  by  being  empirical  but  was  soon  supplemented  from 
analytical  sources.  This  paper  gives  a  survey  of  the  state  of  science  in  the  field. 

It  must  be  admitted  that  the  Importance  of  post-buckling  is  past  its  prime  of 
life.  As  previously  stated,  many  modern  aircraft  structures  are  so  heavily  loaded, 
or  need  so  much  stiffness,  that  their  compactness  yields  high  buckling  stresses. 
Smoothness  of  ’surface  required  for  aerodynamic  per Tormanco  of  high-speed  aircraft  can 
be  another  reason  for  the  expurgation  of  buckled  elements.  However,  some  bright  hopes 
for  the  participation  of  post-buckling  behaviour  in  the  supersonic  age  are  coming 
from  the  complications  to  structural  engineering  caused  by  kinetic  heating  during 
transient  conditions.  Temporarily  high  compressive  stresses  are  being  generated  in 
wing  skins,  which,  added  to  the  aerodynamic  loads,  can  bring  the  skin  into  the  post- 
buckling  range.  Since  this  situation  is  present  only  during  a  small  percentage  of 
the  total  flight  time  (the  heating-up  period),  drag  increase  by  surface  waviness  is 
possibly  not  prohibitive,  privided  this  waviness  does  not  give  rise  to  violent  vibra¬ 
tions.  Aerodynamic  is ts  will  have  to  supply  the  structural  engineer  with  data  on  the 
amount  of  surface  waviness  which  can  be  tolerated. 
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2.  VARIETIES  OF  POST- BUCKLING  BEHAVIOUR 

Buckling  can  occur  only  in  structural  elements  having  at  least  one  dimension 
which  is  much  smaller  than  its  other  dimensions,  as,  for  instance,  in  the  case  of  a 
column,  plate  or  shell  where  the  thickness  is  small  compared  to  the  length  or  width. 
These  structures  are  liable  to  buckling  only  when  their  external  load  is  such  that 
the  internal  forces  act  in  the  direction  of  the  large  dimensions,  for  example,  a 
normal  force  in  a  column,  or  normal  or  shear  forces  in  plates  and  shells.  This  type 
of  load  de'forms  the  structure  in  the  direction  of  the  large  dimensions  and  no  dis¬ 
placement  occurs  in  the  direction  of  the  small  dimension.  In  other  words  the  type  of 
load  which  can  bring  about  buckling  does  not  cause  deflections  before  buckling  starts. 
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At  the  buckling  load  the  state  of  strain  Is  no  longer  unwsbigiously  determined. 
Besides  the  pre-buckling  defornAtion  another  deformation  may  occur,  the  buckling  mode. 
The  characteristic  difference  between  this  deformation  and  the  pro-buckling  deforma¬ 
tion  Is  that  displacements  In  the  direction  of  the  small  size  occur,  e.g.  the  structure 
deflects. 

According  to  the  classic  theory  of  buckling  the  magnitude  of  this  deflection  is 
not  unansbigioasly  determined.  Any  deflection,  provided  It  Is  of  the  first  order  of 
magnitude.  Is  In  equilibrium  with  the  buckling  load.  This  moans  that  the  tangent  to 
the  load-vs-deflection  curve  Is  horizontal  (Pig. 1).  Post-buckling  behaviour  now 
relates  to  the  character  of  the  load-deflection  curve  at  finite  buckling  deflections. 

In  Figure  1  the  various  possibilities  of  behaviour  have  been  schematically  represented. 
Curve  a  applies  to  columns,  where  the  load  is  constant  over  a  wide  range  of  deflec¬ 
tions  provided  no  plastic  deformation  occurs  nt  the  point  of  maxinur  stress.  Due  to 
the  deflection,  the  ends  of  the  column  are  approaching  each  other.  This  end  dis¬ 
placement  Is  a  quadratic  In  the  deflection.  Figure  2a  represents  the  load-vs- 
shortenlng  of  the  column.  With  Increasing  compression  of  the  column  the  end  load 
remains  constant.  As  soon  as  plastic  deformation  starts  the  load  required  for 
equilibrium  decreases,  and  the  column  collapses  (dotted  line  in  Figures  1  and  2). 

This  last  type  of  behaviour  has  something  in  common  with  type  b,  which  applies  to 
shells,  where  with  increasing  deflection  the  load  required  for  equilibrium  decreases. 

It  has  appeared  from  studies  on  the  stability  of  equilibrium  at  buckling  load  that 
the  classic  theory  is  wrong  in  its  statement  that  the  structure  at  the  buckling  load 
is  in  equilibrium  when  the  deflections  are  infinitely  small.  II. L.  Cox2  and  von  Kdrman 
and  Tsien3  discovered  that  some  non-linear  effect,  overlooked  so  far,  was  of  vital 
importance.  The  tangent  to  the  load-deflection  curve  is  not  horizontal";  therefore 
as  soon  as  deflections  occur  under  a  load  equal  to  the  buckling  load  Pt),  the  positive 
difference  between  Ph  and  the  equilibrium  load  accelerates  buckling,  resulting  in 
explosive  collapse.  Only  after  attaining  very  large  deflections  does  the  equilibrium 
load  pass  through  a  minimum.  The  highly  unstable  character  of  the  structure  at  the  load 
Ph  means  that  its  post-buckling  behaviour  is  structurally  useless. 

The  last  possibility  represented  by  curve  c  shows  that  the  equilibrium  load  increases 
with  increasing  deflection.  Since  the  load  increment  as  well  as  the  edge  displacement 
increment  are  proportional  to  the  square  of  the  deflection,  the  load-v3-**dge  displace¬ 
ment  curve  shows  a  positive  slope  at  P^.  This  case  applies  to  flat  plates  supported 
at  their  edges  and  to  slightly  curved  panels. 

So  far  it  has  been  assumed  that  the  structural  members  under  consideration  are 
perfect,  which  means  that  they  do  not  def1 oct  at  loads  below  P^.  When  the  structure  is 
initially  imperfect  -  for  example,  the  column  not  absolutely  straight,  the  plate  not 
absolutely  plane  and  the  shell  showing  analogous  deviations,  it  delfects  before  the 
loa  P^  is  reached.  The  behaviour  of  the  three  types  of  structures  is  illustrated  by 
Figures  3  and  4.  There  is  no  longer  a  critical  load  at  which  buckling  starts. 
for  the  perfect  structures  a  and  b  are  replaced  by  a  load  which  is  close  to  P^  in  the 
case  of  the  column  but  much  smaller  than  in  case  b.  Case  c  on  the  contrary  shows 
a  gradual  increase  of  loan  with  increasing  deflection  or  edge  displacement.  The 
transition  point  between  pre-buck] ing  and  post  buckling  behaviour  vanishes.  Moreover, 
at  loads  far  beyond  P.0  the  curve  approaches  the  curve:  for  the  perfect  structure. 
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Slneo  isper feet  ions  cannot  be  completely  avoided  in  actual  structures  the  liability 
of  type  b  to  ltperfcctions  is  a  warning  to  keep  actual  loads  far  below  P^.  This  puts 
the,  use  of  structures  in  the  buckled  state  still  further  out  of  consideration. 

However,  with  type  c,  imperfections  do  not  alter  the  conclusion  that  these  structures 
are  useful  in  the  post-buckled  range. 

Further  attention  will  be  restricted  to  the  study  of  structure  type  c. 


3.  THE  PHYSICAL  CONCEPT  OF  POST  BUCKLING  BEHAVIOUR 

Consider  a  perfect  structure,  such  ns  a  flat  plate  loaded  up  to  the  buckling  load 
P^.  Then  the  stresser  in  the  plate  are  and  Its  edge  displacements  in  the  plane  of 
the  plate  are  u^.  Now  let  the  edge  displacements  increase  proportionally,  yielding 
u  =  Au^.  Then  the  plate  remaining  flat  has  the  stresses  o  =  Ao^  and  the  required 
edge  load  is  P  =  AP^. 

This  state  of  stress  f.nd  strain  In  unstable.  Allowing  the  plate  to  deflect 
decreases  the  strain  energy  and  therefore  the  undeflected  cannot  be  maintained  and 
the  plate  buckles  ’.n  such  a  pattern  that  the  equilibrium  in  stable.  The  stability 
criteria  for  prescribed  edge  displacement  is  that  the  strain  energy  is  a  minimum; 
any  state  of  strain  adjacent  to  the  stable  one  yields  larger  strain  energy.  Hence 
the  problem  when  investigating  post- buckling  behaviour  is  to  find  for  given  edge 
displacements  the  buckle  pattern  with  the  least  strain  energy.  Denoting  the  strain 
energy  as  A,  this  condition  yields  SA  =  0  for  any  variation  of  the  displacements 
which  is  continuous  and  complies  with  the  kinematic  boundary  conditions.  This  varia¬ 
tional  problem  yields  equations  which  are  identical  with  the  equations  of  equilibrium 
for  the  deformed  element  of  the  plate. 

The  strain  energy  of  the  buckled  plate  consists  for  one  part  of  ‘extenslonal’ 
energy  -  the  work  done  by  the  membrane  stresses  over  the  extensions  of  the  middle 
surface  of  the  plate,  -  and  for  another  part  of  ‘inextensional'  energy  -  the  work  done 
by  the  bending  moments  ovr-  the  curvatures  of  the  deflected  plate. 

In  the  undeflected  state  the  strain  energy  is  completely  extenslonal.  Therefore 
since  the  total  strain  energy  in  the  buckled  state  is  less,  buckling  is  the  way  in 
which  the  plate  gets  rid  of  much  of  its  extenslonal  energy.  This  is  the  general  trend 
governing  the  buckle  patten  preferred  by  the  plate;  at  the  expense  of  some  bending 
energy  the  extenslonal  energy  is  kept  down.  The  more  the  plate  proceeds  in  the  post- 
buckling  range  the  more  the  ratio  of  extenslonal  energy  to  inextensional  energy 
decreases.  The  wave  pattern  tends  towards  a  developable  surface5. 

The  suppression  of  extenslonal  energy  means  that  the  membrane  stresses  are  kept 
down. 

The  statements  are  illustrated  in  Figure  5  by  the  case  of  a  plate  strip  (width  b) 
simply  supported  at  its  longitudinal  edges  and  longitudinally  compressed.  With 
increasing  overall  strain  e  the  relative  magnitude  of  the  membrane  stresses  decreases, 
the  centre  part  of  the  plate  where  the  membrane  stresses  are  small  becoming  wider. 
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This  e^ans  that  the  tuckllng  pattern  flattens  core  and  Bore  with  increasing  strain 
ratio  e/e^.  This  is  evident  free  Figure  6.  Shortening  of  a  longitudinal  fibre  of 
the  plate  is  accomplished  without  excessive  compression  of  its  centre  line  Just  by 
following  a  curved  course.  This  shortening  effect  is  a  consequence  of  the  finiteness 
of  the  deflections  and  is  a  quadratic  in  the  amplitude  of  the  waves.  Since  the 
longitudinal  edges  of  the  plate  are  supported  they  cannot  escape  compression  by 
deflecting  and  they  have  to  carry  large  compressive  stresses. 

Due  to  the  continuity  of  the  deflection  this  edge  effect  extends  over  some  width 
of  plate.  Here  the  membrane  stresses  are  larger  than  tv  buckling  stress.  This 
explains  why  P  exceeds  P^.  A  further  conclusion  is  tha.  che  distribution  of  P  over 
the  plate  width  depends  on  the  degree  of  buckling:  therefore  “.tt*  lysis  of  post- 
buckling  behaviour,  cannot  start  from  given  edge  stresses  but  must  start  from  *lven 
edge  disDlacements  or  from  a  given  edge  load. 

The  favourable  behaviour  of  plates  after  buckling  proves  to  be  a  consequence  of 
the  fact  that  the  edges  are  supported.  If  the  longitudinal  edges  were  free  to  move 
they  would  deflect  as  well  as  the  centre  part  and  the  load  carried  after  buckling 

would  be  constant,  just  as  in  the  case  of  a  column. 

When  the  plate  is  distorted  in  shear  the  same  general  trend  applies.  The  shear 
angle  y  formed  by  the  edges  of  a  rectangular  plate  is  composed  of  two  parts  of 
different  origin  (Fig. 7).  For  one  part  it  is  the  shear  of  the  middle  surface  of  the 

plate  due  to  its  membrane  stress.  For  the  other  part  it  follows  from  the  oblique 

wave  pattern.  By  folding  a  rectangular  plate  in  oblique  waves  the  length  in  >.he 
direction  normal  to  the  wave  crest  is  decreased;  this  makes  the  length  of  the 
diagonals  of  the  rectangle  unequal  and  the  rectangle  changes  into  a  parallelogram. 

Figure  7  does  not  account  for  edge  support  by  which  the  waves  are  being  suppressed  in 

the  edge  region.  Figure  8  shows  how  the  plate  manages  to  make  the  transition  between 
its  straight  edges  and  the  washboard  pattern  of  its  centre  part.  Whereas  in  Figure  7 

the  wave  formation  is  such  that  the  plate  is  still  a  developable  surface  this  is  no 

longer  the  case  when  the  edges  are  kept  straight.  The  main  part  of  the  plate  is  a 
developable  surface,  but  the  edge  region  is  not,  but  near  the  edges  the  membrane 
stresses  vary  rapidly  yielding  a  double  curved  surface.  This  is  evident  from  the 
fact  that  plate  sections  parallel  to  the  sides  of  the  plate  differ  in  developed 
length.  In  tne  case  where  the  plate  is  supported  by  completely  stiff  edge  members 
the  projected  length  of  the  plate  section  does  not  change.  Hence  the  true  length  of 
a  buckled  plate  strip  is  longer  than  that  of  the  edge  strip,  its  extension  being 
p. jportional  to  the  square  of  the  amplitude.  Therefore  the  buckle  pattern  induces 
tensile  stresses  in  the  plate  parallel  to  the  edges,  these  stresses  falling  off 
towards  the  edges. 

If  the  plate  were  to  remain  flat  it  would  sustain  at  the  given  strain  y  the  very 
large  membrane  stress  G>',  yielding  large  extensional  energy.  In  changing  over  to  a 
buckled  shape  it  does  not  completely  get  rid  of  membrane  stresses,  but  it  brings  them 
down  to  a  lower  level  with  correspondingly  smaller  extensional  energy. 

Another  characteristic  feature  of  the  tendency  to  keep  extensional  energy’  down  can 
be  noticed,  namely,  the  constant  level  of  membrane  stresses  over  a  large  part  of  the 
plate.  Since  the  strain  energy  has  something  to  do  with  the  square  of  the  stresses 
and  Ja2 dx  >  /cr|ydx  the  strain  enprgy  is  beirg  reduced  when  cr  keeps  close  to  the 
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average  stress.  Hero  agnin  the  edge  support  is  the  basic  condition  by  which  the 
plate  can  sustain  loads  in  excess  of  l*b.  Without  edge-support  the  panel  could  be 
sheared  according  to  Figure  7  by  oblique  waves.  With  vanishing  bending  stiffness  the 
si  nr  deformation  could  be  obtained  without  any  shear  force.  In  the  case  of  Figure  8 
the  strain  ratio  y/yb  is  large  and  the  plate  has  got  far  into  the  post-buckled  rnnge. 
When  the  ratio  is  close  to  unity  the  buckle  pattern  is  quite  different  and  tends  to 
the  buckle  pattern  at  Pb.  Then  the  surface  has  double  curvature  everywhere  and  the 
developahle  central  part  vunisiies. 

Returning  now  to  general  considerations  some  further  statements  can  be  made. 

Since  at  given  edge  displacements  the  strain  energy  is  smaller  in  the  buckled  state 
than  in  the  flat  state  the  external  load  required  for  this  edge  displacement  is 
actually  smaller  than  in  the  absence  of  buckling.  It  follows  that  Lne  post-buckling 
stiffness,  being  the  ratio  between  load  and  edge  displacement,  i',  smaller  than  the 
pro-buckling  stiffness.  This  is  demonstrated  by  Figure  S.  The  ratio  of  the  incre¬ 
ments  of  load  and  edge  displacement  defines  the  stiffness  of  the  buckled  structure 
with  respect  to  incremental  deformation.  It  is  still  smaller  than  the  post  buckling 
stiffness.  This  incremental  stiffness  is  important  when  considering  the  stability 
of  the  complete  structure,  namely,  general  stability  and  aero-elastic  stability. 

When  loading  a  plate  in  the  post-buckling  it  may  happen  that  the  wave  pattern 
changes  with  explosive  violence.  Obviously  the  preceding  pattern  contained  more 
energy  than  the  final  one,  the  excess  energy  causing  the  bang.  This  seems  to  be  in 
contrast  to  the  statement  that  the  way'  in  which  a  structure  distorts  is  governed  by 
least  energy.  The  following  explanation  of  this  phenomenon  can  be  given. 

With  increasing  edge  displacements  the  buckle  pattern  changes  continuously.  It 
starts  off  with  a  pattern  corresponding  to  the  buckling  mode  and  the  shape  of  the 
wave  modifies  gradually.  The  tendency  toward?  reduction  of  extensional  energy  at  the 
expense  of  more  bending  energy  involves  a  tendency  towards  reduction  of  wave  length. 

The  wavelength  however  cannot  change  freely  due  to  the  edge  conditions  imposed.  For 
instance,  with  a  rectangular  plate  of  aspect  ratio  4  loaded  in  compression  the 
buckling  mode  comprises  4  half-waves.  Thus  number  of  waves  cannot  change  continuosly, 
the  only  choice  being  an  integral  number  of  half-waves.  Therefore,  while  increasing 
the  edge  displacements  the  plnte  sticks  to  4  half-waves,  even  when  at  a  certain  strain 
ratio  5  half-waves  would  yield  less  strain  energy.  The  panel  will  jump  into  the 
pattern  with  lower  energy  level  only  when  its  equilibrium  gets  unstable.  This  will 
not  necessarily  occur  at  the  intersection  of  the  strain  energy  curves  of  Figure  10. 

The  transition  point  depends  on  the  power  of  the  disturbance  which  forces  the  plate 
out  of  its  stable  n-wave  configuration  through  intermediate  configurations  into  the 
stable  configuration  with  more  waves.  Such  disturbances  are  present  in  the  form  of 
imperfections  of  the  plate.  Dependent  on  the  magnitude  of  the  imperfection  the  jump- 
over  occurs  at  a  smaller  or  larger  edge  displacement.  When  the  imperfection  is  very 
weak  it  may  be  that  the  transition  is  retarded  so  much  that  the  (n+l)-wave  pattern  is 
passed  over  and  the  wave  number  changes  from  n  immediately  into  n+2  or  more. 

Summarizing: 

(a)  Edge  support  is  the  basis  of  plate  performance  beyond  buckling. 
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(b)  Analysis  of  pQ«t*buckl lng  behaviour  must  atari  from  given  edge  displacements  nnd 
not  from  given  edge  stresses. 

(c)  The  equilibrium  equations  art*  non-linear  In  the  deflection. 

(d)  For  given  edge  displacements  the  wave  pattern  Is  such  that  strain  energy  is  a 
minimum. 

(e)  The  plate  disposes  of  much  of  Its  strain  energy  by  buckling. 

(f)  Extensional  energy  is  kept  down  at  the  expense  of  inextens tonal  energy  the  more 
the  plate  proceeds  into  the  post-huckl ing  range.  This  reduces  the  membrane 
stresses. 

(K)  The  buckle  pattern  changes  with  increasing  edge  displacements,  tending  to  a  wide 
centre  part  with  constant  small  membrane  stress  and  developable  surface,  to 
narrow  edge  regions  with  large  membrane  stress  and  double-curved  surface, 
and  to  n  short  wavelength. 

(h)  With  finite  plate  length  explosive  changes  of  wave  pattern  towards  increasing 
wave  number  can  occur. 

(i)  The  stiffness  reduces  the  more  the  plate  proceeds  into  the  buckled  range. 


4.  THF.  EQUATIONS  GOVERNING  POST- BUCKLING 
BEHAVIOUR  OF  FLAT  PLATES 


Denoting  by  u,  v,  w  the  displacements  in  the  directions  of  the  axes  of  x,  y,  7. 
respectively,  where  the  axis  of  z  is  normal  to  the  plate,  the  equations  of  equilibrium 
of  a  plate  element,  comprising  the  full  thickness  t,  are 
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Equationa  (la)  and  (lb)  are  quadratic  and  Equation  (1c)  is  cubic  in  w.  When  T 
is  not  constant,  one  coefficient  of  Equation  (lc)  is  a  function  of  x  and  y.  Exact 
solutions  of  these  equations  are  in  general  difficult  to  obtain  and  approximate 
solutions  must  be  established. 


Introducing  the  stress  function  P  defined  by 
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Equations  (la)  and  (lb)  can  be  replaced  by  the  compatibility  condition 
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As  in  the  case  of 


These  equations  are  known  as  the  von  Karman-Tsien- equations 
the  previous  group  they  are  non-linear  and  can  not  be  solved  exactly.  When  edge  dis¬ 
placements  u,  v  are  given,  as  is  usual  in  this  type  of  problem,  the  boundary  conditions 
for  F  must  be  derived  from  the  following  stress-strain  relations: 
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Since  approximate  solutions  can  be  obtained  on  the  basis  of  the  principle  of 
minimum  strain  energy  the  expression  for  strain  energy  per  unit  area  will  be  given. 
It  is* 
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5.  METHODS  FOR  SOLUTION 

Approximate  solutions  always  start  with  some  assumption  about  the  function 
representing  the  deflection  w.  Apart  from  mathematical  reasons  it  is  the  most  spec¬ 
tacular  component  of  the  displacement  and  therefore  the  least  difficult  one  to 
estimate.  The  assumption  of  w  contains  at  least  one  and  preferably  more  undetermined 
parameters  p. 

Cue  method  of  solution  is  along  the  following  lires,  using  the  Equations  (2). 
Substituting  w  into  (2a)  this  equation  is  a  linear  (quation  in  F  and  can  be  solved 
exactly  or  if  necessarily  approximately,  taking  due  account  of  the  boundary  condi¬ 
tions.  Then  F  is  known  as  a  function  of  the  parameters.  Thereupon  w  and  F  are 
substituted  in  (2o).  The  solution  being  an  approximate  one  the  lefthand  side  of  (2b) 
cannot  be  made  to  vanish  b;  choosing  a  finite  number  of  parameters  p  and  some  method 
must  be  devised  for  determining  the  parameters  such  that  (2b)  is  satisfied  in  the 
average. 

For  this  the  Ritz-Galerkin  method  Is  preferable.  It  yields  the  equations 

II  dxdy  ”  0  (4) 

the  number  of  which  is  equal  to  the  number  of  parameters.  If  F  is  the  exact  solution 
of  (2a)  this  method  is  equivalent  to  the  Raleigh-Ritz  method.  It  yields  those  values 
of  the  parameters  p  for  which  the  strain  energy  is  less  than  with  any  other  choice  of 
parameters.  Since  the  solution  is  an  approximate  one  the  energy  is  larger  than  the 
actual  amount  when  the  plate  would  be  free  to  find  the  wave  pattern  corresponding  to 
minimum  energy.  .  Therefore  the  stiffness  obtained  in  this  way  is  larger  than  the 
actual  stiffness. 

It  is  a  pity  that  exact  solutions  for  F  can  be  obtained  only  for  very  simple  functions 
of  w,  and  that  even  then  their  computation  is  tedious8. 

When  F  is  an  approximate  solution  of  (2a)  nothing  can  be  concluded  about  the 
nature  of  the  approximation.  Since  the  compatibility  condition  is  not  satisfied  thp 
strains  are  not  compatible,  whereas  the  Raleigh-Ritz  method  applies  to  compatible 
strains  only.  Approximating  F  by  using  the  minimum  complementary  energy  theorem, 
the  approximation  obtained  for  the  extensional  energy  is  too  small.  This  is  only 
one  part  of  the  ene  gy;  the  other  part  is  being  overestimated  by  using  Equation  (4). 
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So  with  the  mixture  01  the  theorems  of  least  complementary  energy  and  leant  strain 
energy  In  an  tnr  Igoroua  way  no  definite  conclusion  is  possible  on  the  nature  of  the 
approx  ir.nt  ion. 

In  this  respect  a  method  which  is  based  on  one  theorem  only  is  preferable,  the 
theorem  of  minimal  strain  energy.  This  solution  starts  from  Equations  (1)  and  applies 
the  Ritz-Galerkln  method,  cr  it  starts  from  Equation  (3)  and  applies  the  Raleigh- 
Ritz  method.  These  methods  are  mechanical 1>  identical.  Besides  an  assumption  for 
v,  assumptions  have  to  be  made  for  u  and  v.  These  assumptions  include  some  parameters 
p.  Then  the  strains  can  be  computed  and  consequently  from  Equation  (1)  the  'addi¬ 
tional  forces’  qx.  q  ,  q/  or  from  Equation  (3).  the  strain  energy. 

Using  the  Raleigh-Ritz  method,  the  equations  from  which  the  parameters  p  can  be 
solved  are 


JJ^dxdy  =  0  (5) 

Using  the  Ritz-Galerkin  method,  the  equations  for  the  determinations  of  the  parameters 
p  are 


qx,  q^,  qz  are  the  additional  forces  per  unit  area  required  for  equilibrium  at  the 
approximate  displacements  u,  v,  w.  The  condition  expresses  that  the  potential  energy 
of  this  additional  load  system  does  not  change  when  the  parameters  vary.  It  can  be 
shown  that  this  condition  is  equivalent  to  the  minimum  potential  energy  theorem,  which 
forms  the  basis  of  Equation  (5). 

When  u,  v,  w  do  not  satisfy  the  dynamical  boundary  conditions,  equilibrium  requires 
additional  edge  forces.  In  that  case  the  left  hand  side  of  Equation  (6)  should  be 
complemented  by  the  contour  integral  of  the  product  of  this  additional  edge  force  and 
the  derivative  to  p  of  the  corresponding  edge  displacement. 

The  difficulty  with  these  methods  is  that  an  appropriate  guess  of  u  and  v  has  to 
be  made.  Of  course  the  adaptibility  of  the  guess  can  be  improved  by  increasing  the 
number  of  parameters,  but  this  soon  leads  to  a  prohibitive  number  of  non-linear 
simultaneous  algebraic  Equations  of  the  form  (5)  or  (6).  For  this  prectical  reason  u 
and  v  must  be  constructed  by  means  of  the  parameters  adopted  for  w.  The  best  way  to 
do  this  is  to  make  assumptions  on  the  distribution  of  the  strains  e.r,  ev,  y,  from 
which  u  and  v  can  be  derived,  tailing  into  account  the  given  edge  displacements. 

This  procedure  comes  close  to  the  procedure  which  is  upplied  when  (2a)  is  solved 
approximately.  Then  appropriate  stress  functions  F  have  to  be  chosen  and  this  is 
essentially  a  guess  of  stress  distribution  or  strain  distribution.  Then  however  the 
possibility  exists  of  combining  a  certain  number  of  these  stress  functions,  which 
gives  the  approximation  more  adaptability,  provided  the  coefficients  of  the  stress 
functions  become  simple  functions  of  the  parameters  p.  Unfortunately  this  latter 
condition  is  usually  not  satisfied,  which  means  again  that  the  evaluation  of  the 
equations  is  getting  prohibitively  laborious. 
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6.  AVAILABLE  SOLUTION  FOR  BEHAVIOUR  AFTER  BUCKLING 

6.1  Long!  tudinnl  Gxwpresftion  (Effective  Width) 

6.1.1  Flat  plates 

A  longitudinally  compressed  plate  carries  the  load 


P  = 


cr  dy 


fthen  the  edge  strain  (compressive)  is  e  the  stress  in  the  edge  stiffener  is  =  Ee. 
Now  the  ‘effective  width’  bfi  is  defined  by 
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ca.i  be  called  the  post-buckling  effectiveness  of  the  plate. 


Von  K&rman  answered  the  urgent  need  for  knowledge  on  effective  width  by  devising 
and  engineering  guess  based  on  very  much  simplified  considerations.  Although  this 
was  an  ingenious  close  guess,  giving  (Fig. 11) 
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its  reliability  is  not  self-evident. 

The  first  thorough  investigation  is  due  to  Marguerre  and  Trefftz8’9.  They  assumed 
the  wave  pattern 
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where  f  and  L  are  the  parameters.  This  wave  form  comprises  the  buckling  mode,  so  it 
will  yield  a  good  approximation  for  small  strain  ratios  e/e.  .  However,  for  e/e^  =  & 
it  yields  bp/b  =  1/3,  whereas  fo”  physical  reasons  be/b  should  vanish.  This  dis¬ 
crepancy  is  a  consequence  of  expressing  the  wave  form  as  o  function  of  y,  which  does 
not  permit  flattening  at  larger  strain  ratios.  H.L.  Cox10  accounted  for  this  effect, 
assuming  ever  the  centre  part  a  constant  wave  depth  and  sinusoidal  deflection  in  the 
edge  region.  His  parameters  are  the  wave  amplitude  and  the  width  of  tile  edge  region, 
the  wave  length  being  constant.  His  application  of  Raleigh’s  metnod  ir  not  rigorous, 
since  part  of  the  extensional  energy  is  neglected.  In  this  respect  the  energy  Is 
being  underestimated;  however,  the  approximation  of  the  wave  form  has  the  opposite 
effect. 
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The  esoKt  advanced  InvoHtlgatlon  has  been  done  by  W.T.  Kolter11.  Unfortunately, 
this  work  has  only  been  published  In  Dutch,  bo  that  It  Is  little  known.  Kolter 
followed  Cox’  s  nuggesMon  for  the  wave  pattern,  allowing  the  wave  length,  however, 

In  addition  to  the  wave  uepth  and  the  width  of  tne  edge  region,  to  he  a  parameter. 
Also,  he  did  not  neglect  part  of  the  energy. 

Several  assumptions  were  made  for  the  displacement  functions:  for  the  simply 
supported  plate  (4),  one  of  then  having  continuous  9?w/5y*’  at  the  transition  of  the 
centre  part  and  edge  region;  for  the  clamped  plate  (2);  and  for  the  clast tcally 
restrained  plate  (3). 


The  good  agreement  between  the  results  obtained  with  different  wave  forms  indicates 
that  the  approximation  is  satisfactory  and  that  the  results  are  not  sensitive  to 
slight  alterations  In  the  wnve  form  assumed.  Moreover,  the  effective  width  as  a 
function  of  the  strain  ratio  proved  to  be  practically  independent  of  the  boundary 
conditions.  Hence  the  effective  width  for  simply  supported,  elastically  restrained, 
or  clamped  edges  can  be  given  by  the  same  function  of  e/eh,  viz. 


4/  s 
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where  t'b  Is  the  buckling  strain  pertaining  to  the  actual  boundary  conditions. 
Equation  (8)  gives  a  representation  of  this  function  with  no  more  than  1%  error. 
Figure  11  shows  that  the  agreement  between  experimental  results  for  various  edge 
restraints  and  Koiter’s  analytical  results  is  quite  satisfactory.  Since  the  approx¬ 
imation  must  yield  a  too  large  effective  width  it  is  somewhat  surprising  to  find 
that  most  of  the  experiments  fall  above  the  theoretical  curve.  The  most  acceptable 
explanation  will  be  that  in  actual  panels  the  wave  length  reduction  with  increasing 
strain  ratio  will  be  retarded,  whereas  Koiter's  computations  apply  to  the  infinitely 
long  plate,  where  noting  prevents  a  gradual  change  of  wave  length. 


Koiter's  theory  has  the  advantage  that  it  holds  for  any  strain  ratio  between  1  and 
lX.  More  accurate  investigations  h  ive  been  made  by  Hemp1?,  applying  t.o  strain  ratios 
a  little  over  1.  They  are  not  discussed  here,  because  the  inevitable  imperfections 
of  plates  affect  the  behaviour  of  plates  in  this  region  sensibly  (compare  Figure  4), 
reducing  the  practical  importance  of  more  accurate  knowledge. 


During  the  last  few  years  the  Netherlands  Aeronautical  Laboratory  has  proceeded 
with  its  investigations  on  effective  width,  extending  it  to  the  behaviour  of  plates 
in  the  plastic  range  by  experiments  on  24S-T  and  75S-T,  cladded  and  uncladdpd13,  l4. 
The  rnnge  of  strain  ratios  was  1  <  e/t^  <  50  and  Cp/e^  was  in  the  range  1.5  <  Cp/e^ 
<  20,  where  €p  is  the  strain  at  which  the  elastic  limit  of  the  material  is  exceeded. 
Referring  the  effective  width  to  the  edge  stress  =  Eg£,  the  results  suggest  the 
empirical  conclusion  that  Koiter’s  Equation  (8)  may  be  extrapolated  into  the  plastic 
range.  Some  representative  results  are  shown  in  Figure  12. 


Some  series  of  tests  were  also  carried  out  on  2S-&  H  material  so  as  to  investigate 
post-buckling  behaviour  with  materials  having  a  more  gradual  transition  between  the 
elastic  and  plastic  parts  of  the  stress-strain  curve14.  In  this  case  the  effective 
width  in  the  plastic  range  is  greater  than  that  given  by  the  theoretical  curve  for 
elastic  materials.  Some  results  are  shewn  in  Figure  13. 


i:  ■ 
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fi.l.?  iVarrcrer  cylindrical  panels 

Quito  recently  Kotter  investigated  the  Initial  post-buckl  inu  behaviour  of  narrow 
cylindrical  panels  such  as  occur  in  wlnga  and  fuselages  with  longitudinal  stiffeners15. 
All  constraining  effects  of  the  stiffencis  other  than  the  complete  suppression  of 
radial  deflection  were  neglected. 

The  tangent  to  the  post-buckling  curve  at  buckling  load  has  been  established.  It 
is  given  by 

P  -  Pt, 

=  p  +  w  (9) 


where  8  -  [l2(l  -  v2)]'*  b(Rt)‘^ 


<fi(8)  =  1  -  29r'  -  z'Yl  ^(n2  +  1)_2  +  (n2  +  -  1  -  <94]*: 

n=  i 

Pb  =  Ebt(l  +  <9Vb  for  8  $  1 
and  is  shown  in  Figure  14. 

In  wing  panels  the  values  of  8  are  usually  considerably  smaller  than  0. 5  and  ns 
appears  from  Figure  14  the  post-buckling  behaviour  is  very  similar  to  that  for  flat 
plates,  where  8  =  0. 

In  fuselage  panels  the  values  of  8  range  from  approximately  0.4  to  about  2,  so 
that  usually  the  initial  post-buckling  behavioir  differs  considerably  from  that  of  a 
flat  plate. 

For  8  <  0.3  dP/de  is  practically  equal  to  that  for  <9=0;  curvature  does  not 
affect  the  post-buckled  behaviour.  For  8  >  0. 5  dP/de  starts  to  decrease  rapidly  with 
increasing  8,  passing  at  8  =  0.65  through  zero  and  becoming  negative. 

The  horizontal  tangent  represents  the  limiting  case  of  a  stable  post-buckled 
behaviour  for  a  prescribed  load.  At  <9  =  0.77  the  tangent  is  vertical  and  at  8  -  1 
the  behaviour  is  very  close  to  that  of  a  full  cylinder. 

These  figures  apply  to  the  panel  the  edges  of  which  are  simply  supported.  It  may 
‘e  expected  that  initial  post-buckling  behaviour  with  clamped  or  elastically  restraint 
edges  will  be  stable  up  to  higher  values  of  8,  but  the  general  trend  will  be  the 
same. 

The  effect  of  initial  sinusoidal  imperfections  at  small  deflections  has  been 
investigated.  Then  the  lead  up  to  which  the  behaviour  is  stable  is  lowered 
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considerably  when  8  >  0.65  (compare  Figure  4b).  This  'buckling  load’  depends  on  the 
ratio  of  the  amplitude  a  of  the  Initial  tares  to  plate  thickness.  Figure  15  shows  this 
for  two  values  of  0. 

As  for  thj  more  advanced  post-buckling  stages  it  seems  reasonable  to  expect  that 
it  approaches  that,  of  a  flat  plate  of  equal  width. 

Considering  a  free  panel  simply  supported  at  its  edges,  the  Initially  curved  slate 
tends  In  Its  cent  re  part  to  become  a  devolt,  able  surface,  as  with  the  initially  Hat 
plate.  Supposing  that  the  curved  plate  when  buckled  has  the  same  shape  as  the 
initially  flat  plate,  toe  strain  energy  per  unit  length  in  the  first  case  is 
H  Fbt  (£[j  is  the  buckling  strain  of  the  flat  plate)  greater  than  that  of  the 

flat  plate,  where  it  is  n  Ebteb*[-*  +  |  (e/£b)3/j]  in  the  initial  post-buckling  stage 
and  in  the  advanced  stage  even  more  than  that.  Hence  when  e/eb  »  1  the  strain 
energies  in  these  two  cases  aie  approaching  one  another,  yielding  equal  effective 
widths. 

With  a  shell  running  continuously  over  regularly  spaced  stringers  the  individual 
panels  are  not  free  to  assume  the  wave  pattern  of  the  flat  plate,  since  this  would 
imply  a  discontinuity  of  the  slopes  at  the  stringer,  the  discontinuity  being 
h/R  =  2v82e^11 2 .  The  condition  of  continuity  means  that  the  adjacent  panels  are 
restraining  each  other,  the  effect  of  which  is  an  increase  of  stiffness. 

Starting  from  the  buckle  pattern  of  the  single  panel,  the  edges  have  to  be 
rotated  through  the  constant  angle  vd2e^l,!.  This  angle  is  of  the  same  order  of 
magnitude  as  the  slope  of  the  wave  pattern  with  the  free  panel,  when  e/eb  is  small. 
However,  with  a  large  strain  ratio  e/eb  the  slope  with  the  free  panel  is  about  1.7e1/?; 
then  the  correction  required  for  continuity  is  relatively  small  and  so  is  the 
additional  energy  originating  from  this  edge  restraint. 

The  conclusion  is  that  the  load  sustained  by  a  curved  panel  is  for  the  same  edge 
strain  somewhat  greater  than  the  load  carried  by  a  flat  panel  of  the  same  width  and 
tends  to  become  equal  at  large  strain  ratios. 

Figure  16  shows  schematically  the  load-strain  curve  for  plates  which  are  unstable 
in  the  initial  post-huckl Ing  range.  Since  imperfections  cut  off  the  peak  of  the 
load-strain  curve,  near  the  buckling  load  the  curve  for  the  actual  plate  will  be  still 
closer  to  that  for  the  flat  plate.  So  for  practical  purposes  it  is  justified  to 
identify  post-buckling  behaviour  of  narrow  (6  4  1)  cylindrical  panels  and  flat  plates 
in  spite  of  the  initially  unstable  post-buckling  behaviour  when  8  >  0.65. 

6.2  Flnt  Plates  Under  Shear  Load 

As  stated  before,  the  first  analytical  investigation  on  post-buckling  behaviour  of 
plates  loaded  in  shear  was  made  by  H.  Wagner1.  Its  object  is  the  ‘complete  tension 
field’,  where  the  plate  is  supposed  to  be  unable  to  carry  compressive  membrane 
stresses.  This  means  that  yb  =  0,  or  more  precisely  that  the  strain  ratio  (to  which 
this  investigation  applies)  y/yb  =  So  it  represents  the  limiting  case. 

Wagner’s  work  has  been  of  much  importance  in  the  design  of  spar  webs,  etc.,  but  It 
yields  conservative  design  since  the  performance  of  the  plate  is  underestimated.  The 
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analytical  Investigation  of  the  'Incomplete  tension  field’,  where  the  plate  can 
sustain  compressive  as  well  as  tensile  stresses,  starts  with  a  paper  by  Kroran  and 
MarRuerre*  applying  to  the  siaply  supported  plate  strip. 

They  assumed  the  wave  form 


w 


f  sin  v  - 

» 


sin 


v(%  -  my) 


L 


•j)~ 


(10) 


where  f,  L  and  m  are  the  parameters. 

This  wave  pattern  is  an  approximation  of  the  huckle  mode.  Its  nodal  lines  are 
straight,  whereas  the  tangent  of  the  exact  nodal  line  is  perpendicular  to  the  edge. 

The  buckling  strain  yb  obtained  fnm  this  approximate  wave  form  is  5.7%  too  large. 

The  data  obtained  for  post-buckling  behaviour  are  good  at  strain  ratios  not  too 
much  greater  than  1.  At  large  strain  ratios  the  wave  pattern  flattens  In  the  centre 
part,  which  phenomenon  is  not  represented  in  the  pattern  (10).  So  it  is  evident 
that  Reference  9  does  not  yield  Wagner’s  results  In  the  limiting  cace  of  infinite 
strain  ratio.  This  first  attempt  was  improved  by  Koiter56,  assuming,  as  in  the  cose 
of  the  effective  width  problem,  constant  wave  amplitude  in  the  centre  part  (develop¬ 
able  surface)  and  double  curvature  in  the  edge  regions.  This  adds  a  fourth  parameter 
to  those  of  Equation  (10),  namely  the  width  of  the  centre  part. 

Several  assumptions  were  made  for  the  wave  pattern  in  the  edge  regions:  for  the 
simply  supported  plate  (4),  one  of  them  having  continuous  curvature  at  the  transition 
of  centre  part  and  edge  region,  another  one  with  nodal  lines  ending  perpendicularly 
to  the  edges;  and  for  the  clamped  plate  (2). 

This  theory  yields  for  7/>b  =  00  results  which  are  in  agreement  with  Wagner’s 
complete  tension  field  theory.  The  differences  between  the  results  obtained  from 
the  various  wave  patterns  were  not  very  important.  One  of  them  applying  to  the  simply 
supported  plate  was  chosen  for  evaluation  of  Koiter’s  equations,  which  has  been 
established  in  References  17  and  18. 

Reference  17  contains  diagrams  from  which  the  relation  between  the  overall  strain 
ratios  e2/eb’  ”X//fcb  0/1  ^  ^e  averaSe  stress  ratios  a  J Eeb>  aj r/ Et|  can  be 

read.  The  indices  1  and  2  apply  to  the  longitudinal  and  the  lateral  direction 

respectively.  Each  diagram  applies  to  one  of  the  r/E£b  ratics  0.5,  10 .  45.  50; 

the  coordinates  are  ^/e^,  e 2/eb);  the  diagrams  give  curves  of  constant  cr 

cr2/Eeb  and  r/Gy.  Thus  at  given  e  j/e^,  e?/eb  the  corresponding  stress  ratios  a^Ee^, 

o  f Eeb  and  the  corresponding  shear  rigidity  reduction  r/Gy  are  found  by  interpolation 

The  field  of  application  of  these  data  is  not  limited  to  webs  loaded  in  shear. 

Since  e(/eb  and  e?/eb  are  arbitrary  they  apply  as  well  to  the  case  of  stiffened 
plates  loaded  in  shear  and  in  compression  or  tension  longitudinally  and  laterally. 

They  make  it  possible  to  establish  for  simply  supported  plates: 

(a)  Prom  given  edge  displacements  the  edge  loads 
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(b)  Proa  given  edge  load3  the  edge  displacements 

(c)  P'  on  a  mixture  of  three  given  edge  dinplnceinciits  or  loads  the  other  displace¬ 
ments  and  loads 

(d)  The  stiffness  of  the  plate 

(e)  The  stiffness  with  respect  to  load  increments,  for  instance  the  change  of 
shear  stiffness  resulting  from  a  change  of  longitudinal  compression. 

These  data  are  strictly  valid  only  when  the  stresses  are  below  the  proportional 
limit.  The  stresses  can  be  derived  from  the  wave  form.  This  has  been  evaluated  in 
Reference  18,  which  establishes  the  greatest  effective  strain  ep  according  to  the 
Huber-von  Mlses-Hencky  criterion,  again  In  diagrams  applying  to  constant  r/Eeb,  where 
the  coordinates  are  e  /eb,  e2/eb  and  curves  of  constant  eg/eb  are  given. 

The  former  data  apply  to  simply  supported  edges.  The  evaluation  of  plates  with 
clamped  edges  has  not  been  done  systematically.  It  has  been  found  that,  contrary  to 
the  effective  width  problem,  post-buckling  behaviour  of  plates  loaded  in  shear  is 
not  merely  dependent  on  y/y^  but  on  the  boundary  conditions  as  well.  For  the  same 
strain  ratios  y/y^  the  clamped  plate  is  closer  to  the  complete  tension  field  than  is 
the  simply  supported  plate.  This  difference  with  the  effective  width  problem  can  be 
physically  understood  as  follows. 

At  large  strain  ratios  the  strain  energy  of  the  compressed  plate  is  mainly  confined 
to  the  edge  regions  (large  membrane  stresses,  large  curvatures),  whereas  in  the  tension 
field  the  centra  part  as  well  contains  much  strain  energy  since  the  normal  and  shear 
stresses  in  the  centre  part  and  near  the  edges  are  equally  important. 

The  wave  form  flattens  at  some  distance  from  the  edge,  which  means  that  the  effect 
of  boundary  conditions  has  damped  out  when  getting  at  the  centre  part.  Therefore  the 
strain  energy  of  the  centre  part  is  almost  independent  of  the  boundary  conditions. 

On  the  contrary  the  buckling  strain  depends  on  the  boundary  conditions.  Hence  the 
total  strain  energy  cannot  be  merely  a  function  of  the  strain  ratio.  At  very  large 
strain  ratios  the  edge  regions  are  narrow,  so  that  the  total  strain  energy  is  almost 
equal  to  the  energy  of  the  centre  part.  Then  the  strain  energy  would  be  almost 
independent  of  the  boundary  conditions.  This  means  that  for  large  strains  the  post- 
buckling  behaviour  of  clamped  plates  is  almost  the  same  ns  that  of  the  simply 
supported  plate  at  the  same  strain. 

An  experimental  check  of  theory  was  not  possible  with  the  tests  by  Lahde  and 
Wagner10,  since  they  investigated  the  plate  with  clamped  edges.  Tests  on  spar  webs 
are  usually  confined  to  shear  load  and  slight  additional  bending,  and  usually  the 
support  of  the  plate  by  stiffeners  provides  some  restraint.  So  the  Netherland 
Aeronautical  Institute  undertook  a  test  program  on  plates,  loaded  in  compression  and 
shear,  without  restraint  at  the  sup sorting  members.  The  results  have  been  communi¬ 
cated20’21.  The  result  showed  consiierebl p  scatter,  more  in  particular  at  the 
smaller  strain  ratios.  Obviously  the  imperfections  of  the  plate  have  much  effect  on 
post-buckling  behaviour  and  this  scatter  prevented  an  accurate  comparison  of  theory 
and  tests.  However,  the  tests  provide  sufficient  evidence  chat  Koiter' s  theory  gives 
nn  adequate  picture  of  post-buckling  behaviour  over  the  whole  range  of  strain  ratios 
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and  combinations  of  normal  and  shear  loads,  more  in  particular  for  panels  having 
large  aspect  ratio.  In  view  of  scatter  in  post-buckling  behaviour  of  actual  plates 
there  seems  little  point  in  further  refinement  of  analysis.  So  far  Koiter’s  theory 
is  the  only  one  with  this  wide  range.  Several  authors  made  valuable  contributions, 
like  Denke22'2’,  Levy7*’  25  and  Bergmann26,  applying  to  the  square  plate  or  the 
rectangular  plate  with  length  to  width  ratio  not  Inrger  than  2.  However,  the 
validity  of  those  theories  is  restricted  to  small  and  moderate  strain  rntios  and  to 
panels  loaded  in  shear  only. 

The  practical  importance  of  a  satisfactory  theory  after  evaluation  is  that  it 
provides  knowledge  on  stiffness  of  buckled  plates.  From  this  knowledge  the  stress- 
strain  relations  of  a  complete  structure,  composed  of  plates  and  other  elements  like 
stringers  and  frames,  can  be  determined.  A'art  from  that  a  satisfactory  theory 
enables  Jr  to  establish  the  edge  displacements  at  which  permanent  wrinkling  starts. 
However,  theory  is  unable  to  predict  the  edge  displacements  at  which  failure  occurs. 
When  the  structure  fails  in  rivet  joints  before  permanent  wrinkling  starts,  the 
failure  could  be  predicted  by  analysis.  But  usually  rivet  joints  will  be  made 
stronger  and  their  failure  will  be  preceded  by  plastic  deformation  in  the  plate.  In 
view  of  the  complexity  of  the  problem  there  arc  no  hopes  for  analytical  investiga¬ 
tions  on  post-buckling  behaviour  in  the  plastic  range.  Another  possible  cause  of 
failure  is  stiffener  failure  due  to  the  combined  action  of  normal  load  and  flexure 
and  torsion  Imposed  by  plate  wrinkles.  This  phenomenon  too  is  beyond  the  possibili¬ 
ties  of  analysis.  Therefore  knowledge  of  ultimate  lead  of  stiffened  plates  In  the 
post-buckled  range  must  be  acquired  empirically.  In  this  respect  much  has  been  done 
at  N.A.C.A.  by  Kuhn,  Peterson,  Levin  and  other.  Kuhn  has  summarized  this  work  in 
his  recently  published  text  book27.  These  data  enable  us  to  predict  ultimate  loads 
with  no  more  than  10%  error. 

6.3  Flnt  Plntes  Under  Tbeminl  Stresses 

In  supersonic  flight  the  temperature  of  the  structure  rises  due  to  kinetic 
heating.  When  flight  has  been  maintained  at  constant  speed  during  some  length  of 
time  the  temperature  differences  In  the  structure  are  small  and  so  are  the  thermal 
stresses.  However  in  transient  conditions  heat  flow  produces  temperature  differ¬ 
ences  and  consequently  thermal  stresses,  which  affect  both  buckling  and  post-buckling 
behaviour.  The  skin  is  heated  by  the  air,  but  at  the  plate  edges  part  of  the  heat 
input  is  transferred  to  the  interior  structure,  such  as  webs  of  multi-web  wings. 
Therefore  the  edges  have  a  lower  temperature  than  the  centre  part.  Since  the  longi¬ 
tudinal  strain  of  the  plate  strip  is  constant  throughout  the  width  of  the  plate  and 
the  resultant  force  should  vanish,  compressive  stresses  occur  in  the  centre  part  and 
tensile  stresses  occur  near  the  edges.  To  this  stress  distribution  is  added  the 
stress  due  to  the  mechanical  loads  on  the  aircraft. 

Since  we  are  interested  in  post-buckling  behaviour,  the  part  of  the  skin  which 
should  be  examined  is  the  side  which  is  loaded  in  compression.  Here  the  skin  plate 
presents  before  buckling  a  stress  distribution  of  the  type  represented  in  Figure  17. 
It  has  been  shown20  that  a  good  approximatikn  of  the  buckling  mode  is  given  by  the 
buckling  mode  for  uniformly  stressed  plates.  This  is  a  guide  for  post-buckling  in¬ 
vestigations.  So  the  buckling  pattern  appropriate  to  non-thermal  conditions,  assumed 
in  Reference  11,  can  be  applied  to  thermal  conditions. 
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In 


In  Reference  20  post. -buck  11  ns  behaviour  has  been  investigated  for  plates  loaded 
compression  and  by  the  following  temperature  distribution: 


T  =  Tu  v  T  Mn  77  y/b 


(11) 


or  by  a  temperature  distribution  in  which  T  is  constant  ove'  part  of  the  plate  width 
in  the  centre  and  falls  off  townrds  the  edges  according  to  a  power  n  (Fig.  18). 
Conclusions  from  this  investigation  can  be  given  in  a  very  simple  form,  since  the 
load-strain  curve  appears  to  have  the  same  shape  as  that  for  the  non-thermal  condi¬ 
tion.  Denoting  by  P  and  e  the  load  and  compressive  edge  strain  for  the  non-thermal 
condition,  the  relation  between  P  and  c  in  the  thermal  condition  follows  from: 
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The  edge  stress  is  ap  =  E(e  +  aTe). 


(12) 


So  the  P-e  curve  is  as  shown  in  Figure  19;  it  follows  from  the  P-£  curve  simply 
by  sliding  the  axes  over  a  distance,  which  depends  on  the  temperature  distribution. 
The  buckling  load  decreased  by  thermal  stresses.  However  the  load  increment  in  the 
post-buckled  range  is  independent  of  the  temperature,  it  depends  only  on  the  incre¬ 
ment  of  edge  strain  beyond  the  critical  strain.  This  statement  is  incorrect  in  so 
far  It  does  not  consider  the  effect  of  temperature  on  the  modulus  of  elasticity.  It 
must  be  complemented  by  the  Instruction  that  the  P-e  curve  refers  to  the  modulus  of 
elasticity  for  the  temperature  level  under  consideration. 


The  integral  over  T  in  Equation  (12)  shows,  that  the  buckling  mode  across  the 
width  is  proportional  to  sin  77y/b.  When  allowing  for  flattening  of  the  centre  part, 
as  applied  in  Reference  11,  the  P-e  curve  is  only  slightly  different  for  e/e^  >  4. 
With  the  temperature  distribution  (11)  the  temperature  effect  is  rather  pronounced; 
nevertheless  the  difference  between  P  computed  from  (12)  and  P  computed  from  the 
flattened  wave  shape  was  well  within  1%.  Equation  02) ,  therefore,  can  be  considered 
to  hold  for  the  whole  range  of  strain  ratios. 


The  simplicity  of  the  conclusion  that  the  relation  between  load  and  strain  incre¬ 
ment  is  almost  unaffected  by  thermal  conditions  suggests  that  some  physical  evidence 
must  exist.  The  following  explanation  can  be  given. 


Consider  two  equal  plates  1  and  2  differing  only  in  temperature  distribution  -  the 
temperature  of  plate  1  being  constant,  plate  2  having  varying  temperature  -  longitudinally 
compressed  up  to  their  respective  buckling  loads. 

At  the  buckling  load  Pb  the  stresses  are  crb.  Since  the  strain  energy  is  a  minimum 
for  any  variation  by  of  the  strain  y,  which  Is  accompanied  by  the  longitudinal  edge 
displacement  L&£, 
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/  X(ab&>')dv  =  PbLi>f  (13) 

for  both  plates  1  and  2,  where  <?b  and  Pb  are  different.  Now  consider  the  plates  both 
to  be  compressed  logni tod lnol ly  the  some  amount  Ac  beyond  their  compressive  strain 
eb  and  assuae  that  the  wnvo  pattern  of  plate  2  Is  the  same  as  that  of  plate  1.  e.g, 
the  assumption  Is  that  the  wave  pattern  Is  affected  negligibly  by  the  difference  ot 
thermal  conditions.  Then  the  stress  increments  Ao  following  from  Ac  are  equal  for 
plates  1  and  2.  The  load  increment  Is  AP.  Since  plate  1  is  in  equilibrium  In  Its 
deflected  shape 

/  I[(ob,  *  Aa.tSyldv  ,  (Pbl  ♦  Api)L6£  (!4a) 

for  any  compatible  system  of  strain  variations  Sy.  The  approximation  for  the  lond 
upon  plate  2  Is  found  from  the  condition  that  the  work  done  by  P  is  equal  to  the 
Increase  of  strain  energy  when  e  increases  with  Sc.  The  strain  variations  between 
the  states  of  compression  c  and  e  +  Sc  are  Sry.  Hence 

/2[(obj  +  Aap&yjdv  =  (Pb?  +  Ap?)L5e  (14b) 

Substituting  (13)  into  (14a)  and  (14b),  and  remembering  that  A oy  =  A o?  it  follows 
that  Apj  =  A  P;  for  equal  compressive  strain  increment  At  beyond  buckling  strain. 
This  conclusion  has  been  based  upon  the  assumption  that  the  wave  pattern  of  the  non- 
thermal  case  Is  a  good  approximation  for  the  wave  pattern  in  the  thermal  case.  This 
restricts  the  validity  of  the  statement  that  A  P/A  6  is  unaffected  by  thermal  condi¬ 
tions,  but  the  impression  is  that  these  restrictions  is  not  important. 

Our  reasoning  was  given  with  respect  to  the  longitudinally  compressed  plate.  How¬ 
ever,  any  other  edge  displacement  could  have  been  taken.  Heive  we  can  extend  the 
validity  of  the  statement  to  any  case  of  thermal  load  and  mechanical  load  in  shear 
and  compression,  provided  the  wave  pattern  is  affected  only  slightly  by  thermal 
stresses. 

This  widens  the  statement  to:  the  relation  between  load  increment  P-Pb  and  overall 
strain  increment  e-eb  is  not  influenced  by  thermal  stresses  When  the  post-buckling 
behaviour  Is  known  for  non-thermal  conditions,  the  problem  of  post-buckling  for 
thermal  conditions  is  solved  as  soon  as  the  buckling  load  Pb  has  been  established. 

For  the  time  being  the  only  available  literature  on  post-buckling  behaviour  of 
thermally  loaded  plates  Is  contained  in  Reference  7,  The  case  considered  there  was 
devised  so  that  experimental  check  of  the  analysis  was  possible. 

Since  the  temperature  distribution  is  far  from  actual  conditions  with  kinetic 
heating  of  aircraft  the  results  need  not  be  discussed  in  detail.  It  suffices  to 
say  that  the  theoretical  and  the  measured  wave  patterns  agreed  satisfactorily. 
Furthermore,  Reference  7  establishes  an  interesting  method  for  computing  the  correc¬ 
tion  of  the  deflections  following  from  initial  Imperfections,  which  are  of  the  sR-ne 
form  as  the  wave  pattern. 
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Schematic  representation  of  equilibrium  load  versus  buckling  deflection 
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Schematic  representation  of  equilibrium  load  versus  edge  displacement 
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Fig. 5  Relative  decrease  of  membrane  stresses  with  increasing  edge  displacement 
beyond  buckling.  Longitudinal  compressive  strain  e  of  a  long  rectangular 

plate  of  width  b 


Pig. 6  Shortening  of  a  plate  due  to  buckling 


Shear  angle  of  a  rectangular  plate  following  from  oblique  waves 
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8  Wave  pattern  for  large  strain  ratio  yf 7^ 
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Stiffness  of  buckled  panels 


epresentation  of  post-buckling  behaviour  in  advanced  stage  for 
narrow  cylindrical  panels 


